Magnetization plateaus of the Shastry-Sutherland model for SrCu 2 (B03)2: 
Spin-density wave, supersolid and bound states 



O 

o 
o 

(N 

> 

o 
O 



c3 



I 

O 

o 



(N 
> 
O 

o 

O 
O 

o 



: 

o . 
o ■ 



X 



Tsutomu MomoiD 

Lyman Laboratory of Physics, Harvard University, Cambridge, MA 02138 

Keisuke Totsuka 

Department of Physics, Kyushu University, Hakozaki, Higashi-ku, Fukuoka-shi, 812-8581 Japan 

( ) 

We study the Heisenberg antiferromagnet on the Shastry-Sutherland lattice under magnetic fields, 
to clarify the magnetic properties of SrCu2(B03)2- Treating magnetic excitations promoted by 
the field as Bose particles and using strong coupling expansion, we derive an effective Hamiltonian 
for the effective magnetic particles. Anisotropic repulsive interactions between effective particles 
induce 'insulating' states with a stripe SDW structure at magnetization m/m sa t = 1/3 and a 
checkerboard structure at 1/2, and thereby form magnetization plateaus. Supersolid phases appear 
around insulating SDW phases by changing the magnetic field. Nature of these supersolid phases is 
discussed in detail. We also demonstrate how the geometry of the Shastry-Sutherland lattice affects 
dynamical properties of magnetic excitations significantly and makes a novel type of quintuplet 
(S = 2) boundstates condense for very small magnetization. 



I. INTRODUCTION 

Since plateau structures were observed in the magne- 
tization process of a series of quasi one-dimensional Ni- 
compoundsEl, magnetization plateaus have been attract- 
ing extensive interests. The appearance of plateaus in 
magnetization curves was explained as metal-insulator 
transitions, of magnetic excitations driven by a mag- 
netic fielcu; Magnetic excitations crystallize and form 
SDW in the plateau states, and they are itinerant in 
the non-plateau states. Recently it was discussed that 
this phenomenon is not limited to the one-dimensional 
systems, but more general, and. occurs in two- and three- 
dimensional systems as wellHTJ. Before the recent stud- 
ies, magnetization plateaus were already known to ap- 
pear at m/m sa t = 1/3 in the antiferromagnetic com- 
pounds on the triangular lattice, CgEuD, CsCuChjQ, 
and RbFe(Mo04)2P. Theoretically, plateaus were seen 
at m/m sat = 1/3 in the pHeisenberg antiferromagnet 
on the triangular latticeBilil, and also at m/m sa t = 
1/2 in the multiple-spin exchange model with four-spin 
interactionsMa. The l/3pplateau comes from appearance 
of a collinC|M|4tMrf stataalij and the 1 /2-plateau from the 
uuud stateOa. These magnetization plateaus can also 
be regarded as superflukLinsulator transitions of Hipped- 
spin degree of freedomcffl. It may be worth mentioning 
that there is also another trial to realize magnetization 
plateau in two-dimensional systems as gappedjspin liquid 
states analogous to the FQHE wave functionsEj. 

Recently a quasi two- 

dimensional compound SrCu2 (1303)2 is attracting exten- 
sive interests because it shows magnetization plateaus 
and peculiar dynamical properties. The two-dimensional 
lattice structure of Cu 2+ ions in i%iGu2 (603)2 is the 
so-called Shastry-Sutherland latticeEJiia, which is shown 
in Fig. 0. Susceptibility and specific heat data show 



that interactions are antiferromagnetic, the spin exci- 
tation has a gap above the ground state, and the spin 
anisotropy is weaktJ. This material seemstj to be well 
described with the S — 1/2 HeisepJperg antiferromagnet 
on the Shastry-Sutherland latticeEH. (Hereafter we call 
this model simply as Shastry-Sutherland model.) The 
ground state of the Shastry-Sutherland model is exactly 
a direct product of local dimer singlets on bonds J for 
the region J' I J < 0.68 (Refs. HIJ) and there is a fi- 
nite gap above the ground state. Susceptibility and spe- 
cific heat estimated from this model .with J'/ J = 0.68 
fit well with the experimental results&3. (Recently the 
value has been modified to J'/J = 0.635 by taking 
into account the three-dimensional coupling of Shastry- 
Sutherland layers £3) In ref. |l5|, Kageyama et al. reported 
two plateaus at m/m sat = 1/8 and 1/4 in the magnetiza- 
tion curve of SrCu2 (603)2- Theoretically, we studied the 
magnetization process of the Shastry-Sutherland model, 
treating a dimer triplet as a particle, and thereby pre- 
dicted a novel iproad plateau at m/m sat = 1/3 in the 
previous reports. It was argued that the appearance is 
due to superfluid-insulator transition of the excitations. 
Quite recently, the above 1 /3-plateau was experimentally 
observed in magnetization measurements up to a strong 
field of 57[T]lii As was predicted in ref. [5], this plateau 
is the broadest one ever found in this material. This 
seems to support the correctness of our argument based 
on the particle picture. Xhis_material also shows pecu- 
liar dynamical propertiesEj~E3, e.g. one-magnon excita- 
tion is almost dispersioplpss, but two-magnon excitations 
have strong dispersior£3. The aim of this paper is to 
present the details of our analyses and results reported 
briefly in ref. |^, and to proceed further thereby giving 
remarkable consequences of the correlating hopping of 
the effective Hamiltonian. This correlated hopping can 
also explain the peculiar dynamical behaviors observed 
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FIG. 1. Shastry-Sutherland lattice. Solid (dotted) lines de- 
note bonds with strong exchange J (weak one J'). 



in experiments^. 

The ground state of the Shastry-Sutherland model at 
zero magnetic field was studied Jar a varying ratio J'/ J 



exact diagonalization 
It was foundu that 



by the mean field approximatior 
methocO, and series expansion 
the ground state is the exact dimer state for J'/ J < 0.69, 
a gapped plaquette singlet state for 0.69 < J'/ J < 0.86, 
and Neel-ordered state for 0.86-,-j^-, J'/ J. The parame- 
ters of SrCu 2 (B0 3 ) 2 estimatedNlia as J' /J = 0.68 (or 
0.635) suggest that the spin state of the real material be- 
longs to the dimer phase, but it is very close to the phase 
boundary with the plaquette singlet phase. Consistency 
between theoretical and experimaalal results on magne- 
tization plateau at m/m sat = 1/3EE3 f>*j£. dynamical be- 
havior in inelastic neutron scatteringEJnil also supports 
that the real material is in the dimer phase. 

In this paper, we study the 5=1/2 Heisenberg antifer- 
romagnet on the Shastry-Sutherland lattice and discuss 
the magnetic properties under a magnetic field. We ana- 
lyze this model using strong-coupling expansion. In sec- 
tion [n| we derive an effective Hamiltonian for the dimer- 
triplet excitations. Virtual triplet excitations yield vari- 
ous effective repulsive interactions, which are responsible 
for the plateaus. Although the usual single-particle hop- 
ping is completely missing from the resulting Hamilto- 
nian, correlated hopping processes are contained instead. 

In section Eljj, we investigate the magnetization pro- 
cess applying the classical approximation to the effective 
(pseudo-spin) Hamiltonian and show that plateaus ap- 
pear at m/m sa t = 1/3 and 1/2. Both of the plateau 
states are Mott (SDW) insulators of magnetic excita- 
tions. Spatially anisotropic interactions perturbatively 
generated stabilize several SDW structures. Near the 
plateau states, there are supersolid phases, in whichjSDW 
long-range order (LRO) and superfluid coexist£3~E3. 
Field (B) versus coupling (J'/ J) phase diagram is also 
presented. 

Although the inter-dimer coupling J' /J is not small, 
the special geometry of the Shastry-Sutherland lattice 
strongly suppresses bare one-particle hopping and hence 
the correlated hopping becomes important. This fact 
leads to remarkable consequences on the motion of mag- 
netic excitations. In section IV, we investigate correlated 



hopping more closely and demonstrate how it favors the 
formation of a bound pair of dimer triplets which repel 
each other at the level of the bare Hamiltonian. This 
bound state has a relatively large dispersion and may be 
an elementary particle at very low magnetization, instead 
of the single dimer triplet. The non-plateau state would 
be superfluid of these bound states at least for very small 
magnetization. Above a certain threshold value of mag- 
netization, individual dimer triplets become elementary 
particles and the non-plateau state is characterized by 
superfluidity of single dimer triplets. 

In section we also discuss spin excitation just below 
the saturation field. It is found that lowest energy states 
exist on a close curve instead of a point in the momentum 
space. 

Critical phenomena of the plateau transition are dis- 
cussed in section VI and are argued to be in the same 
universality class as the superfluid transition of the in- 
teracting boson system in the dilute limit. 



II. EFFECTIVE HAMILTONIAN 

In this section we derive an effective Hamiltonian for 
the magnetic excitations under a strong enough magnetic 
field. We begin with the J' = limit. In this limit, the 
lowest triplet excitation over the (dimer) singlet ground 
state is apparently obtained by promoting one of the 
dimer singlets to a triplet. Although the dimer prod- 
uct remains to be the exact ground state even for non- 
zero J', the above completely localized triplet does not; 
perturbation J' 'broadens' the triplet by exciting nearby 
singlets. Unlike the ground state which is perfectly free 
from quantum fluctuation, excited states strongly suffer 
from it. This is one of the most important features of the 
model. As a result, (physical) triplets can interact with 
each other with the help of virtual triplets created by per- 
turbation and the effective Hamiltonian for the physical 
triplet degrees of freedom should contain effective inter- 
actions. 

The most systematic way to take into account such ulc.- 
tual processes would be the strong-coupling expansionBfl. 
We start from the limit J' = and treat the interaction 
with coupling J' by perturbation. In the absence of the 
external field, the spin states of a single isolated dimer 
(J) bond consists of a singlet and triplet separated by a 
gap J. When the field is increased until the lowest triplet 
with S z = 1 intersects the singlet, we may keep only two 
states-the singlet and the lowest triplet-as the physical 
degrees of freedom as far as the low-energy sector is con- 
sidered. 

We carry out the degenerate perturbation for such 
low-energy degrees of freedom. Considering the triplet 
(S = 1) state with S z — 1 as a particle (a hard-core 
boson) and the dimer singlet (S = 0) as a vacancy, we 
derive an effective Hamiltonian for the magnetic parti- 
cles. (Rest of spin states, i.e. triplets with S z — and 
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— 1, are included into the intermediate virtual states of 
the perturbation.) The perturbational expansion is per- 
formed up to the 3rd order in J' /J for degenerate states 
with a constant number of dimer triplet excitations with 
S z = 1. The final form of the effective Hamiltonian is as 
follows 



H — Hq + H i + H 2 + i?3 , 
H = {J-B)Y,nu 



J' x , 

H x = — ^riiUj, 

{id) 
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+ 2 ( & l+e2+eA + h.c.)n i+e2 n l+( 
+2(bl_ e2 _ el bi + h.c.)ni-e2rii-ei 
+ 2 (C-e2+eA + h.c.)ni-e 2 n i+ ei 
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where indices i and j run over an effective square lattice 
of dimer bonds (both horizontal and vertical), and A (B) 
sublattice contains horizontal (vertical) ones. The oper- 
ator b\ (pi) creates (annihilates) magnetic particle with 
spin S z = 1 at bond i, and n% — b\bi. The unit vectors 
ei and e 2 are shown in Fig. [l]. The interactions summed 
up on B sublattices are obtained by replacing ei and e 2 
in those on A sublattices. The abbreviations NN and 
NNN denote pairs of nearest neighbor and next nearest 
neighbor sites. 

The Hamiltonian derived above does not have the bare 
one-particle hopping terms like b\bj + bjbi, which means 
a single dimer triplet excitation is dispersionless at this 
order. This is in a striking contrast with other spin gap 
systems e.g. the 2-leg ladder. The energy gap of one 
dimer triplet excitation is evaluated as AEi (B) = J — 
B-J' 2 /J~J' 3 /(2J 2 ). 

On the other hand, the effective Hamiltonian con- 
tains many correlated- hopping processes, where an ef- 
fective hopping of a particle is mediated by another 
one. Roughly speaking, these are closely related to 2- 
particle Green functions of the triplet bosons. This cor- 
related hopping is important for dynamics of the Shastry- 
Sutherland model and it leads to many interesting con- 
clusions; bound states without attraction, supersolid 
with stripe structure, etc. Most terms of higher orders 
concern the correlated hopping. 

Many two-body repulsive interactions are also de- 
rived, whose range and geometry are shown in Fig. 4 of 
Ref. p. The nearest neighbor repulsion Vnn and the next- 
nearest-neighbor one Vnnn are derived as Vnn = J'/2 + 
J' 2 /2J - J /3 /8J 2 and Vnnn = J /3 /4J 2 . (There were a 
typographical error and a mistake in the coefficients of 
the 3rd order terms of Vnn and Vnnn in ref. ||. There 
was also the same mistake in ref. ^7].) The 3rd- neighbor 
repulsion V3 rc i between vertical (horizontal) bonds with 
the distance 2ei (2e 2 ) is V 3rd = J' 2 /2J + 3J' 3 /4J 2 . Thus 
V3 r d is anisotropic and acts only in one direction. The ef- 
fective Hamiltonian (1) correctly reflects the space-group 
symmetry of the original Shastry-Sutherland model. If it 
is considered as a model of interacting hard-core bosons 
on a square lattice consisting of dimer bonds, it still con- 
tains matrix elements which are not invariant under naive 
90° rotation about one site. We show in section III that 



they lead to anisotropic SDW states with a stripe struc- 
ture. Longer-range repulsions between particles can ap- 
pear from higher-order perturbations. Actually, Miya- 
hara and UedaEZI independently took into account such 
long-range interactions in a phenomenological manner, 
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also adding weak one-body hopping by hand. Our find- 
ing is that there is much stronger correlated hopping pro- 
cesses and that triplets are not necessarily localized. 



III. MAGNETIZATION PROCESS 

In order to discuss the ground state properties of the 
effective Hamiltonian ([!]), we first consider only repul- 
sive interactions. From naive consideration of the range 
and geometry of the repulsive interactions (see Fig. 4 in 
Ref. g), we can imagine various insulating density- wave 
states of dimer triplets. Let us treat the repulsive inter- 
actions obtained at different orders of J'/ J, separately. 

(a) The strongest interaction comes from the 1st- 
order term and it is repulsive for adjacent 
triplets. This repulsive interaction chooses SDW 
with 2-(dimer)sublattices checkerboard structure 
at m/m sa t = 1/2, where the unit cell has 4 sites in 
the original Shastry-Sutherland lattice. (See Fig. 
6(b) in Ref. |.) 

(b) The next strong repulsive interactions originate 
from 2nd order terms and they are contained in 
Vnn an d Vnnn ■ Because of anisotropic interaction 
V^rd, a stripe SDW state with a 6-(dimer)sublattice 
structure is stabilized at m/m sa t = 1/3 as shown 
in Fig. [| 

(c) The 3rd-order perturbation generates the weakest 
interactions, which 
favor SDW of 10-(dimer)sublattices checkerboard 
structure at m/m sa t = 1/5. 

These SDW states do not necessarily appear if we con- 
sider all interactions and correlated hopping terms to- 
gether. Whether SDW insulator is realized or not is de- 
termined by competition between repulsive interactions 
and correlated hopping. 

To consider both effects of correlated hoppings and in- 
teractions together, we study the effective Hamiltonian in 
the classical limit. To this end, we map the hard-core bo- 
son system to the S — 1/2 quantum (pseudo)spin system 
and then approximate the Pauli matrices by the compo- 
nents of a classical unit vector. Within the mean-field 
approximation, we first search for the ground state tak- 
ing into account two and four (dimer)sublattices with a 
checkerboard structure and also three (dimer)sublattices 
with a stripe structure, because insulating states with 
these configurations are expected from the above con- 
sideration. To take account of larger sublattice struc- 
tures, we next study the ground state of a finite sys- 
tem with a Monte Carlo method, where we gradually 
decrease temperatures to zero. We consider a finite clus- 
ter of 60 dimers and impose a periodic boundary condi- 
tion that matches with all SDW configurations expected 
from the repulsive interactions, i.e. 2-(dimcr)sublatticc 



and 5-(dimer)sublattice checkerboard structures and 3- 
(dimer)sublattice stripe one. 

The evaluated magnetization processes are shown for 
the cases J/J' = 0.45, 0.6, and 0.68 in Figs. |(a), (b), and 
(c), respectively. Plateau structures appear at magneti- 
zation m/m sa t = 1/3 and 1/2, but no plateau appears 
at m/™ sat = 1/5 and 1/4. Near rn/m sat = 1/4, the 
slope of the curve becomes less steep, but not flat. This 
means that the m/m sat = 1/4 state is energetically sta- 
ble, though it does not have spin gap. In the plateau 
phases at m/m sat = 1/3 and 1/2, there arc SDW or- 
ders and no off-diagonal long-range order (ODLRO), that 
is, collinear spin states are realized. In the non-plateau 
phases, spins have ODLRO. Spin configuration of each 
phase is explained in the following subsections. A re- 
mark is in order here about our classical approximation. 
We obtained several plateaus by analyzing the classical 
pseudo-spin Hamiltonian. From this, one may conclude 
that these plateaus are of a classical origin. However, the 
spins approximated by vectors are not the original ones 
but are pseudo spins obtained for quantum objects-spin 
singlet and triplet-and these plateaus are not classical. 

As has been mentioned above, the spin configurations 
obtained in our approximate method concern the pseudo 
spin defined by triplet (S z — +1) and singlet on a dimer 
bond. In order to translate the classical pseudo-spin con- 
figuration into the original S = 1/2 one, it is convenient 
to consider the so-called spin-coherent state which real- 
izes almost 'classical' states using quantum states. It is 
well-known that an S = 1/2 coherent state specified by 
a classical unit vector 



Q = (cos $ sin O, sin $ sin O, cos O) 



(6) 



is given by (aside from a phase factor coming from the 
gauge degrees of freedom) 

|$, 6) = e*i* sin ( |singlct) + e" 2 ** cos (^j (triplet) 

(7) 

Plugging the expressions of | singlet) and | triplet) in terms 
of original S = 1/2, we obtain 



0\ 1 



$,e>=e^*sin^j- 7 =(|U)-|IT)) 
+e-^*cos(|W). 



(8) 



Taking expectation values 
1 



(*,e|5+|$,e> 



2V2 
1 



e i<E> sin© 



($,e|S+|$,e} = +-^e**sin© 
($,91^1$, 6} = ($,e|5||$,e) = i(l + cosO) , (9) 
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(b) J'/J = 0.6 




FIG. 3. Finite-size (N — 24) cluster that is used for the 
exact diagonalization study. 

we can see that in the presence of the superfluid LRO 
(O 7^ 0, 7r) the off-diagonal elements of the original 
5=1/2 spins on a dimer bond align in an antiparal- 
lcl manner. 

In order to check the accuracy of the strong coupling 
expansion and the classical approximation, we also stud- 
ied a finite system of the original Shastry-Suthcrland 
model using exact diagonalization method. We diagonal- 
ized 24-sites system with a periodic boundary condition 
(Fig. ||) that matches with the 12-sublattice struc ture a t 
magnetization above m/m sat = 1/3 (see section HI B| ). 
The results are shown in Fig. |^. Total behavior shows 
good agreement with the results from the strong coupling 
expansion. 

In the following, we discuss the nature of each phase 
and phase diagram. 



A. SDW on plateau states 



1.2 




0.5 1 1.5 2 2.5 3 3.5 4 
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On the plateau states, dimer triplet excitations crystal- 
lize forming SDW long-range orders, and there are spin 
gaps. The plateau states at m/m sat = 1/3 and 1/2 have 
SDW with stripe and checkerboard structures, respec- 
tively, which are shown in Figs. 6(a) and (b) in Ref. [|. 
These configurations are consistent with the insulating 
states naively expected from the range and geometry of 
repulsive interactions. The particles are perfectly closed 
packed at m/m sa _ t = 1/2 and 1/3 avoiding the repulsive 
interactions from 1st- and 2nd-order perturbations, re- 
spectively. (See Fig. [| for the case of m/m &&t — 1/3.) 



B. Supersolid on non-plateau states 



FIG. 2. Magnetization processes of the Shastry-Sutherland 
model with J'/J = 0.45 (a), J'/J = 0.6 (b), and J'/J = 0.68 
(c). Solid line is obtained from the effective Hamiltonian in 
the classical limit and dotted one from the exact diagonaliza- 
tion (ED) of the original Shastry-Sutherland model with the 
size N = 24. The solid line for J'/J = 0.68 (c), taken from 
Fig. 5 of Ref. 0, is also shown for the comparison with ED. 



By applying a stronger magnetic field than the critical 
value, the plateau states continuously change to super- 
solid states in which superfluid components appear and 
coexist with SDW of the plateau states. Since the ap- 
pearance of the superfluid component is accompanied by 
the Goldstone bosons, the plateau gap collapses. 

Above the magnetization m/m sat = 1/3, particles with 
density 1/3 still form SDW with the stripe structure, 




FIG. 4. Stripe SDW structure at m/m sat = 1/3, which 
is stabilized by the anisotropic repulsion Vard- Black bonds 
denote dimer triplet excitations and white ones dimer singlets. 
Arrows show range of repulsive interactions. 



which equals to the one at m/m sa t = 1/3, and the rest of 
particles become superfluid in 'canals' between stripes. 
(See Fig. |(a).) Particles can hop across a line of SDW 
with the help of correlated hopping terms (see Fig. ||) 
and this hopping makes correlations between superfiuids 
in canals. These extra particles can Bose-condense and 
the phases of the superfluid particles align fcrromagneti- 
cally inside an individual canal and antiferromagnetically 
between canals. 

The supersolid state above m/rn sat = 1/2 maintains 
the same SDW as the m/m sa t = 1/2 case and have su- 
perfluid components as well. Phases of superfluid compo- 
nents form a stripe structure and the phase of one stripe 
aligns antiparallelly to those of next ones. (See Fig. ^(b).) 

In both supersolid phases, a single dimer triplet can 
move by itself assisted by SDW LRO and it behaves as 
an elementary particle of the superfluid. The SDW forms 
a global network of crystallized particles and it helps a 
dimer triplet to hop along the network owing to the cor- 
related hopping. We roughly evaluate the hopping of an 
extra particle by treating the stripe SDW as a classical 
background. In Fig. we show the matrix element of 
single-dimer-triplet hopping along a stripe. Note that 
triplet excitations are not confined in a one-dimensional 
space ('canal') between two rows of particles. A dimer 
triplet can hop both parallel and perpendicular to the 
stripe. The superfluid component of the system may be- 
have as an anisotropic two-dimensional Bose gas. The 
hopping matrices inside of a stripe are negative, whereas 
those across a line of triplet excitations are positive. This 
makes the phases of superfluid component ferromagnetic 
inside the stripe and antiferromagnetic between stripes. 
A similar anisotropic superfluid was reported for. the 
Bose-Hubbard model with frustrating interactional. 
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FIG. 5. Spin configurations in supersolid phases slightly 
above m/m sat = 1/3 (a) and 1/2 (b). Black bonds denote 
dimer triplet excitations and white ones with arrows denote 
linear combination of dimer triplet and singlet, which can be 
written as eq. (|^). The SDW components equal to those at 
m/m aa ,t = 1/3 and 1/2, respectively. Superfluid components 
of magnetic excitations exist on white bonds, and arrows on 
white bonds denote phases $ of particles (pseudo spins) . The 
original spin configuration is recovered by eq. (^). 



C. Phase diagram 

We show phase diagram for J' / J vs. B/J at zero 
temperature in Fig. |7[ where the phase boundaries are 
not very accurate. The plateau at m/m sat = 1/2 ap- 
pears only in the region < J'/ J < 0.51, and one at 
m/m sat = 1/3 in < J'/ J < 0.95. For large J'/ J, in- 
sulating phases disappear and become superfluid phases. 
This is because correlated hoppings are dominant in the 
higher order terms and they become efficient for large 
J' /J. The phase diagram may be not qualitatively accu- 
rate for large J'/ J, since our arguments are based on a 
strong coupling expansion. Furthermore there is a possi- 
bility that the elementary particles change to plaquette 
tripletscj for large J'/ J, which is not taken into account 
in our approximation. 

Spin configurations at various values of magnetization 
are summarized as follows: 
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\ / 
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FIG. 6. Hopping processes of a single dimer ex- 
citation (shaded one) near a line of triplet excita- 
tions (black bonds). ti = (J'fn/AJ + 3( J') 3 n/8J 2 , 
ta = ~{J') 2 n/U - 5(J') 3 n/16J 2 , t 3 = -( J') 3 n 2 /8J 2 , and 
U = (J'fn 2 /8J 2 . Here we assume the particle density on 
each black bond to be n. Matrices of hopping processes across 
the line (ti, £4) are all positive and other ones are negative. 
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FIG. 7. Phase diagram at T = with the parameter J'/J 
and the magnetic field B/J. The lower- and upper bold lines 
correspond to the 1st critical field and the saturation field, 
respectively. For the latter, an analytic expression is available 
(see section V). Inserted figures denote density- wave order 
and phases of particles. Black bonds denote localized particles 
(dimer triplets) and white ones vacancies (dimer singlets). 
Arrows on dimer bonds denote phases of particles. 



For < m/m sat < 1/3, 

spin configuration has large sublattice structures. 
The system is highly frustrated and we sometimes 
reached to different ground states with different 
magnetizations in the Monte Cairo method. We 
expect that weak interactions or quantum fluctua- 
tions can drastically change the ground state prop- 
erties. One reason for this ambiguity may be that 
a single dimer triplet is not an elementary parti- 
cle at very low magnetization, but a bound pair of 
triplets is. (See section IV) 



• At m/m sat = 1/3, 

there is only SDW LRO with stripe structure. (Fig. 

§ 

• For 1/3 < m/m sat < m c i with 0.5 < J'/J < 0.8, 
and for 1/3 < m/ra sat < 1/2 with J'/J < 0.5, 

the ground state shows supersolid with stripe struc- 
ture. (Fig. |(a)) 

• At m/m sat = 1/2 for J'/J < 0.5, 

there is only SDW LRO with checkerboard struc- 
ture. (Fig. 6(b) in Ref. 0.) The results of the exact 
diagonalization (Figs. 0(a) and (b)) indicate that 
this phase seems to be more stable than in our ap- 
proximation because of quantum fluctuations. The 
phase boundary at J'/J = 0.51 might move to a 
larger value due to quantum effects. One possible 
reason for this discrepancy is that our classical ap- 
proximation favors wave-like states (e.g. superfluid) 
against the localized state. 



For 1/2 < m/m s&t < m c2 , if J'/J < 0.5, 

fins become supersolid with stripe structure, 
b)) 



(Fig. 



IV. 



For m c i < m/m sa t < 1 with 0.5 < J'/J, and for 
m C 2 < m/m sa _t < 1 with J'/J < 0.5, 
large sublattice structures, e.g. a helical structure, 
appear. We will discuss nearly saturated region in 
Section 0. A new quantum phase may appear in 
this region. 



BOUND STATE OF TWO DIMER TRIPLETS 
NEAR M = 



In this section, we consider a striking effect of corre- 
lated hopping on the dynamics of the Shastry-Sutherland 
model. It is most clearly seen in the low-magnetization 
region where the number of magnetically excited triplets 
is small. 

First we suppose that there are only two excited 
triplets. When they are far aparl from each other, an 
individual triplet can hardly hoplij and it gains little en- 
ergy by moving on a lattice. (This almost localized prop- 
erty was actually observed in inelastic neutron-scattering 
experimentsc3.) On the other hand, when the two are 
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FIG. 8. Four bases |i, 1) (1), \i,2) (2), |i,3) (3), and \i, 
(4) of a bound state of two dimer triplets. 
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9. Chemical unit cell of the Shastry-Sutherland lat- 



adjacent to each other, the situation is completely differ- 
ent. From the effective Hamiltonian, we can easily see 
that correlated-hopp'mg processes make coherent motion 
of two triplets possible, where a pair of triplet dimer ex- 
citations form a bound state with S z — 2. (In the same 
way we can easily derive various bound states with S = 0, 
1, 2 at zero magnetic fieldcJ. Here we only discuss the 
state with S z = 2, which becomes dominant under the 
magnetic field.) 

Using the effective Hamiltonian, we can exactly show 
bound states. Because of the correlating hoppings, one 
triplet excitation exists necessarily in nearest neighbor or 
next nearest neighbor sites of the other. A little calcu- 
lation shows that the hopping processes are decomposed 
into the center-of-mass motion and the relative motion, 
and that the latter is closed within the four states shown 
in Fig. H. Hence we can write the pair excited states as 



^ exp(iP • r 2 ){ Cl (P)|i, 1) + c 2 (P)|i, 2) + c 3 (P)|z, 3) 

+c 4 (P)|M>}, (io) 

where the two-dimensional momenta P = (p x ,p y ) are de- 
fined with respect to the chemical unit cell of the Shastry- 
Sutherland lattice and unit vectors are defined in Fig. |j. 
(In this sense, the readers should not be confused our P 
with that used in refs. |34| and |2(].) The energy spectra 
w 3 (P) (j = I, 2, 3, 4) of them are computed by diagonal- 
izing the following hopping matrix: 



/ 2Vq + Vnnn 



where 



2Vq + Vnn 
— JNNe lPx 



J^rd 

2Vb + Vnn 



\ 

-J N Ne ' lPx 

2Vo + Vnnn j 
(11) 



Vq = J — B — 



J 



2J 2 



Jnn = 


J' 2 


5J' 3 




4J " 


1 9 

16J 2 




Vnn = 


f 

y + 


J' 2 
23 


J' 3 
8J 2 


Vnnn = 


j' 3 

4J 2 "' 






^3rd = 


J' 2 


3J' 3 

\ 7T. 




4J " 


8J 2 





(12) 



In addition, there is another type of correlated motion ob- 
tained from the above one by reversing the space about 
y axis. Roughly speaking, this corresponds to the inter- 
change of dimers A and B, and its spectra are given by 
t^ji^PxiPy)- These branches together with a dispersion- 
less band of a pair isolated triplets give the entire spectra 
of two-triplet sector. 

Since these bound states can move because of corre- 
lated hopping process, energy spectra of these states are 
dispersive. We show only the lowest branch in Fig. [!(]. 
The lowest energy is given at P = (jr, it) as 

AE 2 (B) = 2 V + - {- J 3rd + Vnn + Vnnn 



-\J (^3rd — Vnn + Vnnn)^ 



16 I 2 



(13) 



On the other hand, two independent dimer triplet ex- 
citations (scattering state) have the dispersionless en- 
ergy 2AEi(B) —. 2Vp (small dispersion appears at 6th 
order and higheiOo). Expanding the right-hand side 
of eq.jl^) in J' /J, we can readily verify that there is a 
gain in kinetic energy of - J( J'/ J) 3 /4 from the 2-particle 
threshold. Actually, the lowest energy of the bound state 
is smaller than that of two independent dimer triplets 
state for any J'/J as shown in Fig. [fl]. For example, for 
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FIG. 10. Energy spectrum of the effective Hamiltonian up 
to 3rd order at J'/J = 0.63. A two-triplet bound state con- 
sists of four branches and we show only the lowest one (see the 
upper band), which is apparently lower than the two-particle 
threshold. The spectrum of a single dimer triplet excitation 
(Vo) is also shown (lower band. Note that it is dispersionless 
at this order.) 



J' / J = 0.68 and B = 0, the dispersion of a bound state 
takes the minimal value 0.442 J at P — (tt,tt), whereas 
two independent dimer triplets have the energy 0.761J 
in total. As is easily seen, two triplets combined to form 
the S = 2 bound states actually feel repulsive interactions 
between each other (i.e. Vnn,Vnnn > 0). Since strong 
repulsion Vnn acts for a pair on adjacent bonds, one may 
naively expect that such bound motions are not energet- 
ically favorable. However, we can take an optimal linear 
combination of the four relative configurations \i, a) with 
a = 1, 2, 3, 4, so that the bound state may avoid the effect 
of Vnn (note that only |i, 2) and \i,3) feel the repulsion 
Vnn) while gaining the kinetic energy by the coherent 
motion. For large Vnn we can easily verify that the ef- 
fect of Vnn is canceled in Ai?2- The energy gain due to 
the motion is larger than the cost from repulsion, and 
hence relatively stable bound states are formed. 

Next, we apply a magnetic field in the z-dircction. 
Because S z — 2 bound states have lower energy than 
two unbound triplets, the bottom of the bound states 
first touches the ground state when the field is increased 
and one-triplet excitation still has a finite energy gap at 
the critical magnetic field. By increasing the magnetic 
field more than the critical field, a macroscopic num- 
ber of bound states condense instead of the single dimer 
triplets. Thus the non-plateau state at very low magneti- 
zation may be a superfluid of bound states and is different 
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FIG. 11. Energy of two dimer triplets for various J'/J. 
Dashed line denotes the energy of two non-interacting dimer 
triplets (2Vb) and solid line the energy of the bound state of 
two dimer triplets with the momentum (n,n). 



from the almost localized dimer triplet state discussed 
by Miyahara and UedatLS. The critical magnetic field of 
magnetization process, where the magnetization starts to 
increase from zero, corresponds to half of the energy gap 
A£?2(0). If we regard, the slow increase of magnetization 
below Hi = 22.5[T]E_3 as the consequence of bound state 
condensation, we can estimati^ihe gap Ai? 2 (0) as 51.0[K] 
from magnetization processEjO. On the other hand the 
one-triplet energy ,-gap A£'i(0) is estimated as 34. 7 [K] 
frorrusjisccptibilityE_3, inelastic neutron scattering^!], and 
ESRcil. If we set the parameters as J — 81.4[K] and 
J' = 53.5 [K], estimates of the energy gaps, A2?i(0) and 
A£^(0), coincide with the experimental results. 

We expect that these bound states are destroyed at 
higher magnetization. Indeed, one-triplet excitations can 
mov e arou nd above m/m sa t = 1/3, as discussed in sec- 
tion III B , and they can gain more kinetic energy than 
bound states. There must be a transition of elemen- 
tary particles from bound states for very low magneti- 
zation to one-triplet states for high magnetization. Wc 
can see this transition, if we neglect interactions be- 
tween particles and apply mean-field approximation to 
the correlated hopping term. As particle density is in- 
creased, one-particle hopping process effectively appears 
as (b\bj + H.c.)(n) from the correlated hopping term and 
then triplet particles gain kinetic energy. In this rough 
estimation, a pair of unbound triplets have lower en- 
ergy than the bound state above m/m sa t = 0.253 for 
J'/J = 0.63. We can expect that bound states disin- 
tegrate above a finite value of magnetization and one- 







triplet particles turn to elementary particles. This esti- 
mate of critical magnetization can be highly modified by 
strong correlation and readers should not consider the 
above value seriously. 

Finally we mention about crystallization of bound 
states. One may expect bound states to crystallize at 
low magnetization, but it will not occur. If bound states 
tend to crystallize, bound states loose kinetic energy and 
then they will be unbound by repulsive interaction be- 
tween triplets. 



V. NEAR SATURATION 

There is another region where quantum effects man- 
ifest themselves. In this section, we briefly discuss the 
region just below the saturation field. In this region the 
one-particle excitation can be obtained from the original 
Shastry-Sutherland model without any approximation. 

The single-particle excitation over the fully polarized 
ground state is given by a single flipped spin. The dis- 
persion of this excitation is readily computed by diago- 
nalizing the following hopping matrix 



/-i(J + 4J') 4(1 + 



(1 



-iPy \ 



V 4(1 + e l Py) 



i(J + 4J') 4( 1 + e_lPy ) 2 



4(1 + e'Py) -±(J + 4J') 4(1 + e i *») 



(l + e-**-) -±(J + 4J') / 
(14) 



Because of a dispersionless mode on the closed curve, 
the density of states (DOS) near saturation magneti- 
zation is like ID and different from the usual 2D one. 
Hence the singularity of magnetization curve near satu- 
ration magnetization would be ID-like, i.e., |m c — m\ ~ 
-J\H C — H\, This behavior can be seen in the results of 
the exact diagonalization in Fig. |[ 

The location of the dispersion minimum p min together 
with the corresponding eigenvector determines the spin 
structure at the semiclassical level. Usually spin states 
just below thnxptical field are correctly given by the clas- 
sical modeScll Detailed analyses of the wave functions 
reveal that the fourfold-degenerate classical helical or- 
der, which was pointed out by Shastry and Sutherland!!!!, 
corresponds to the four apexes (±p max , 0), (0, ±p max ) 
(Pmax = 2cos _1 (J'/J)) of the closed curve. The fate 
of the classical helical LRO when the quantum effects 
are fully taken into account is unclear at present. The 
four ground states with helical LRO are connected by 
a gapless line, which means stiffness of helical order is 
vanishing. We hence expect that quantum fluctuations 
destroy the classical LRO. 

On the other hand, for J'/ J > 1 the saturation field 
H c = 4J' becomes J-independent and the band mini- 
Wim shrinks to a single point p = (0,0); the spin struc- 
t ure realized in the vicinity of saturation is the classical 
staggered one where spins connected by J'-bonds align 
ntiferromagnetically in the xy-plane. Correspondingly, 
the transitipn to saturation is the same as that of 2D 
superfluicO. 



In the above equation, momenta (p x ,Py) are defined with 
respect to the chemical unit cell. Reflecting the fact that 
a single unit cell contains four spins, the spectrum con- 
sists of four bands. Note that the four eigenvalues are 
invariant under the point group D 2( j: 

(Px,Py) l-> (Py,Px) and (p x ,p y ) H-> (-Px,Py) • (15) 

In the dimer limit J'/J <C 1, the lower two correspond 
to energy of a singlet particle on dimer bonds, where 
the number two comes from the two mutually orthogonal 
dimer bonds, and the higher two to a triplet (S z — 0) one; 
our approximation in section [n| corresponds to neglecting 
the latter as higher- lying. 

Although the expression of the dispersion is rather 
complicated, we can locate the position of the band min- 
ima in the momentum space, which is relevant in deter- 
mining the structure in the vicinity of saturation. For 
< J'/J < 1, the lowest band takes a minimal value 
— J(l + J' / J) 2 on a closed curve 



cosp x + cos p y 



This implies that magnetization saturates at 



H c = J 1 



j 



(16) 



(17) 



VI. CRITICAL PHENOMENA 

According to an analogy to many-particle theories, a 
plateau state corresponds to a SDW insulating state and 
gapless ones to supersolids. As the plateau states at 
m/m sa t = 0, 1/3 and 1/2 collapse by increasing the 
magnetic field, superfluid components appear, whereas 
the SDW structure still remains in each phase. Con- 
trary to one dimension, superfluid LRO is accompanied 
by the gapless Goldstone mode in our two-dimensional 
case. Hence, collapse of plateau occurs at the same point 
with the onset of superfluidity component. On the other 
hand if we look at the lower phase boundaries of the SDW 
phases, SDW structures change at the transition points 
in our approximation. In this case the phase transitions 
are presumably of first order. 

Now let us consider the case of 2nd order transition. 
We can imagine two different situations. The first one 
is (i) the transitions driven by changing the external 
field while the coupling J'/J is kept fixed. This tran- 
sition is seen in the actual magnetization process. This 
type of field-induced transitions has a close relationship 
to the filling-control insulator-to-superfluid transitions 
in Bose systems, and the methods used there can be 
imported to our case. The second one is (ii) continu- 
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ous quantum phase transitions occurring with fixed val- 
ues of magnetization. As the limit J'/ J f oo is ap- 
proached, the Shastry-Sutherland model reduces to the 
ordinary 5=1/2 square-lattice Heisenberg antiferro- 
magnet, where no magnetization plateau appears. Hence 
magnetization plateaus vanish at some critical values 
of J' I J and are superseded by supersolid (or, super- 
fluid) phases. In the present model, there is no particle- 
hole symmetry around insulating phases apparently. We 
hence conclude that the above two transitions have the 
same universality class. 

First of all, we have to keep in mind that because of 
the special geometry of the Shastry-Sutherland model 
there is no a priori reason for believing that the sys- 
tem is described by the ordinary Bose liquid with a well- 
defined one-particle dispersion e(p) ~ p 2 . Actually, the 
low-order effective Hamiltonian (|l|) lacks the one-particle 
part. However, we have seen in section IV that this leads 
to the formation of dispersive two-triplet bound states 
in the low-field region, and in section IIIB that SDW 
structure makes one particle dispersive in the supersolid 
phase around m/m sat = 1/3 and 1/2. 

Near the phase boundaries, the superfluid amplitude 
is small and we can map the problem onto the effective 
Ginzburg-Landau model described by the superfluid or- 
der parameteiO. (These Bose particles are not necessar- 
ily dimer triplet excitations, but they can be plaquette 
triplet states of two dimer triplets or flipped spin statesd.) 
This enables us to conclude that superfluid-onset transi- 
tion would be described by the (2+l)-dimensional clas- 
sical XY model when the particle-hole symmetry exists 
and-tw-jthe z = 2 mean-field-like transitions when it does 



notBEii In our case, the effective boson model in section 
II does not have particle-hole symmetry and we hence 
conclude that the plateau transition is of the dynamical 
exponent z = 2 and behaves as 

|m - m c | ~ \H - H c \ log(C/\H - H c \) 

(2D system). (18) 

in two dimensions. Here m c and H c denote the critical 
magnetization and field at the plateau transition. In the 
real material, there are weak interactions between two- 
dimensional layers and these interactions will push the 
system above the upper critical dimension, i.e. 

\m - m c | ~ \H - H c \ 

(3D and quasi-2D systems). (19) 

Note that these forms ape quite different from that in 
one-dimensional systems,El \m — m c \ ~ \/\H — H c \. 

Here we also give a comment about the case the kinetic 
term of Bose particles does not behave as e(p) ~ p 2 . 
According to the scaling argument of Ref. [58[ exponents 
of this kind of transitions should satisfy the relation 

zv = 1 . (20) 

If the one-particle dispersion is of the form e(p) ~ p l (I 
denotes an even integer), the standard power-counting 



argument shows that the upper critical dimensions are 
given by d u = I. For d < I, we obtain |m — m c \ ~ 
\H-H c \ d ' 1 . 

For the case that SDW structure changes at phase 
boundary, e.g. the boundary between the plateau state 
at m/m sat = 1/2 and the supersolid one below m/m sat = 
1/2, the system shows a 1st order phase transition in our 
approximation. We don't find any incommensurate phase 
between them. To take into account the possibility of any 
2nd-order phase transition or incommensurate phase, we 
need to consider the effect of quantum fluctuations more 
seriously. 



VII. DISCUSSIONS AND FUTURE PROBLEMS 

To summarize, we studied the magnetic behavior of 
the Shastry-Sutherland model using strong-coupling ex- 
pansion. Magnetic excitations show insulator-supersolid 
transitions at magnetization m/m sa t = 1/3 and 1/2, and 
thereby create magnetization plateaus. Magnetization 
curve obtained neaMp/m sat = 1/3 looks similar to the 
experimental result Ji2l 

At zero magnetization, bound states of triplet exci- 
tations are formed by the correlated hopping process. 
Above the critical field, quintuplet (5 = 2) bound states 
become elementary particles in the ground state and they 
condense. Whereas, for large magnetization, the bound 
states are destroyed in the ground states and triplet exci- 
tations become elementary particles. Triplet excitations 
are essential for the plateau transitions at m/m sa t = 1/3 
and 1/2. In the experiments, it is unclear in which re- 
gion bound states appear as elementary particles. One 
possibility is that they appear in the tail of the magne- 
tization process below m/m sa ± ~ 0.025. In this region, 
slope of the-magnctization curve is different from that of 
rest parts.E§E3 Further detail analysis are needed on the 
nature of quasiparticles. 

In the present analysis, we did not find the plateaus 
at m/m sa t = 1/8 and 1/4. The mechanism of stabiliz- 
ing these plateaus is not yet clear. It may be natural to 
believe that dimer triplet excitations are crystallized by 
longer-range repulsive interaction. (5 = 2 bound states 
cannot crystallize because of a special origin of the bind- 
ing energy as we discussed in section |y|.) We can con- 
sider two origins for the repulsions as follows: 

1) Though we cut the perturbation series at 3rd order, 
the higher-order expansions can produce longer- 
range repulsions between particles. These repul- 
sive interactions may induce crystallization at low 
magnetization. 

2) Longer-range repulsions can come from other anti- 
ferromagnetic spin interactions in the original spin 
model, which have not been accounted for in the 
pure Shastry-Sutherland model. If we treat the an- 
tiferromagnetic interactions between next-neighbor 
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dimer bonds, we can produce crystallization at low 
magnetization. For example, Miiller-Hartmann et 
al. demonstrated the appearance of 1/4-plateau 
considering another spin interaction, which acts be- 
tween nearest-neighbor dimers.c3 

It is unclear which spin interaction is important in the 
real material. We need first-principle calculation to esti- 
mate exchange couplings. We also need to keep in mind 
that the real material is close to the plaquette singlet 
phase. Under the magnetic field, if this phase becomes 
more stable than the dimer singlet state and plaquette 
triplets become elementary particles, plaquette triplets 
may crystallize and hence create magnetization plateaus 
at low magnetization as discussed in refs. j5|]26|. This 
should be considered in a future. 
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APPENDIX A: EFFECTIVE HAMILTONIAN 

In this appendix, we briefly explain how to obtain the 
effective Hamiltonian in the framework of degenerate per- 
turbation. We suppose that the ground states of the 
unperturbed Hamiltonian (Ho) are degenerated and we 
diagonalize these degenerate ground states by degener- 
ate perturbation. Let the operators Po and Qo be the 
projection operators onto the degenerate (unperturbed) 
ground-state sector and its complement, respectively. In 
addition, we define a projection P onto the perturbed 
ground-state sector. 

According to ref. Ell the problem of degenerate per- 
turbation reduces to solving the following problem 



PqHPP \GS. 



= E a K\G.S. 



(Al) 



where the hermitian operator K is defined by 

K = P PPo . (A2) 

However, this form is not so convenient to our pur- 
pose because it does not take the form of the .ordinary 
eigenvalue problem. A trick invented by BlochcS solves 
this difficulty. The key is to introduce a new operator U 
through the following relation 



PP Q = UK . 
The operator U can be expanded as 



u = j2 xnjj(n) . 



(A3) 



(A4) 



where the n-th order coefficients U^ n ' are given by 

U<® = P a 
U (1) = SVP 

U (2) = (SVSV + S 2 VS°V) P Q 
u {3) = (svsvsv + S 3 VS°VS°V + S 2 VSVS°V 
+SVS 2 VS°V + S 2 VS°VSV) Po 

(A5) 

In the above equation, we have used a short-hand nota- 
tion S defined by 



S k = 



-P for k = 

-Qo for k > 1 . 



(E-H ) k 

Then, eq. ( |A1[ ) can be recasted as 
(P HU - E a )K\G.S.; a) = 







(A6) 



(A7) 



This is a usual eigenvalue problem and the matrix we 
have to diagonalize is finally given by 



H c ff = PoHU 

= Eo + XPoVPo + \ 2 PoV 
1 



+A 3 
~PoV 



P V 



Eo — Ho 
1 



E 
QoV- 



Ho 
1 



(Eo — Ho) 



Eo — Ho 
2 Q VP VP 



QoVPo 
QoVPo 



(A8) 



Note that non-hermitian terms appear in general when 
we proceed to terms higher than second order. Reality of 
the eigenvalues is no longer guaranteed. To remedy this 
shortcomings, we have used the average of H e g and H\ s , 
which is now hermitian. 



* On leave of absence from Institute of Physics, University 
of Tsukuba, Tsukuba, Ibaraki 305-8571, Japan. 
Y. Narumi, M. Hagiwara, R. Sato, K. Kindo, H. Nakano, 
and M. Takahashi, Physica B 246-247, 509 (1998); 
Y. Narumi, R. Sato, K. Kindo, and M. Hagiwara, 
J.Mag.Mag.Mat. 177-181, 685 (1998), and references cited 
therein. 

2 K. Totsuka, Phys. Rev. B 57, 3454 (1998). 

3 A. K. Kolezhuk, Phys. Rev. B, 59, 4181 (1999). 



12 



Sakamoto, and K. Kubo, Phys. Rev. B, 59, 



4 T. Momoi, H. 
9491 (1999). 

5 T. Momoi and K. Totsuka, Phys. Rev. B 61, 3231 (2000). 

6 H. Suematsu, K. Ohmatsu, K. Sugiyama, T. Sakakibara, 
M. Motokawa and M. Date, Solid State Commun. 40, 241 
(1981); T. Sakakibara, K. Sugiyama, M. Date and H. Sue- 
matsu, Synthetic Metals, 6, 165 (1983). 

7 H. Nojiri, Y. Tokunaga and M. Motokawa, J. Phys. (Paris) 
49 Suppl. C8, 1459 (1988). 

8 T. Inami, Y. Ajiro, and T. Goto, J. Phys. Soc. Jpn. 65, 
2374 (1996). 

9 H. Nishimori and S. Miyashita, J. Phys. Soc. Jpn. 55, 4448 
(1986). 

10 A.V. Chubukov and D.I. Golosov, J. Phys. Cond. Matter 
3, 69 (1991). 

11 A. E. Jacobs, T. Nikuni and H. Shiba, J. Phys. Soc. Jpn. 
62, 4066 (1993). 

12 K. Kubo and T. Momoi, Z. Phys. B 103, 485 (1997). 

13 Yu.E. Lozovik and O.I. Notych, Solid State Commun. 85, 
873 (1993). 

14 R.W. Smith and D.A. Keszler, J.Solid State Chem. 93, 430 
(1991). 

15 H. Kageyama, K. Onizuka, Y. Ueda, N.V. Mushnikov, T. 
Goto, K. Yoshimura, and K. Kosuge, Phys. Rev. Lett. 82, 
3168 (1999); H. Kageyama, K. Onizuka, T. Yamauchi, Y. 
Ueda, S. Hane, H. Mitamura, T. Goto, K. Yoshimura, and 
K. Kosuge, J. Phys. Soc. Jpn. 68, 1821 (1999). 
S. Miyahara and K. Ueda, Phys. Rev. Lett. 82, 3701 (1999). 
S. Shastry and B. Sutherland, Physica 108B, 1069 (1981). 
S. Miyahara and K. Ueda, unpublished. 
K. Onizuka, H. Kageyama, Y. Narumi, K. Kindo, Y. Ueda, 
and T. Goto, J. Phys. Soc. Jpn. 69, 1016 (2000). 
H. Kageyama, M. Nishi, N. Aso, K. Onizuka, T. Yosihama, 
K. Nukui, K. Kodama, K. Kakurai, and Y. Ueda, Phys. 
Rev. Lett. 84, 5876 (2000). 

H. Nojiri, H. Kageyama, K. Onizuka, Y. Ueda, and M. 
Motokawa, J. Phys. Soc. Jpn. 68, 2906 (1999). 
T. Room, U. Nagel, E. Lippmaa, H. Kageyama, K. 



Onizuka, and Y. Ueda, |cond-mat/9909284 
P. Lemmens, M. Grove, M. Fischer, G. Giinthero dt, V 
N. Kotov, H. Kageyama, K. Onizuka, and Y. Ueda, cond- 



mat/0003094 . 
K. Totsuka, S. Miyahara, and K. Ueda, preprint. 

25 M. Albrecht and B. Mila, Europhys. Lett. 34 145 (1996). 

26 A. Koga and N. Kawakami, Phys. Rev. Lett. 84, 4461 
(2000). 

27 S. Miyahara and K. Ueda, Phys. Rev. B 61, 3417, (2000). 

28 A. F. Andreev and I. M. Lifshitz, Zh. Eksp. Teor. Fiz. 56, 
2057 (1969) [Sov. Phys. JETP 29, 1107 (1969)]. 

29 A.J. Leggett, Phys. Rev. Lett. 25, 1543 (1970). 

30 K.S. Liu and M.E. Fisher, J. Low. Temp. Phys. 10, 655 

(1973) . 

31 M.E. Fisher and D.R. Nelson, Phys. Rev. Lett. 32, 1350 

(1974) . 

32 H. Matsuda and T. Tsuneto, Suppl. Prog. Theor. Phys. 46, 
411 (1980). 

33 E. Roddick and D. Stroud, Phys. Rev. B 48, 16600 (1993); 
A. van Otterlo and K.-H. Wagenblast, Phys. Rev. Lett. 
72, 3598 (1994); 

GG. Batrouni, R.T. Scalettar, G.T. Zimanyi, and 



A.P Kampf, ibid. 74, 2527 (1995). 

34 Z. Weihong, C.J. Hamer, and J. Oitmaa, Phys. Rev. B 60, 
6608 (1999) 

35 E.G. Batyev and L.S. Braginskii, Zh. Eksp. Teor. Fiz. 87, 
1361 (1984) [Sov. Phys. JETP 60, 781 (1984)]. 

36 S. Gluzman, Phys. Rev. B 49, 11962 (1994). 

37 S. Doniach, Phys. Rev. B 24, 5063 (1981). 

38 M. P. A. Fisher, P. B. Weichman, G. Grinstein, and D. S. 
Fisher, Phys. Rev. B 40, 546 (1989). 

39 D. S. Fisher and P. C. Hohenberg, Phys. Rev. B 37, 4936 
(1988); 

V.N. Popov, Functional Integral and Collective Excitations, 
Chap. 6, (Cambridge University Press, 1987). 

40 E. Miiller-Hartmann, R. R. P. Singh, C. Knetter, and G. 
S. Uhrig, Phys. Rev. Lett. 84, 1808 (2000). 

41 T. Kato, Perturbation Theory for Linear Operators, 
(Springer- Verlag, Berlin, 1976). 

42 C. Bloch, Nucl. Phys. 6, 329 (1958). 



13 



